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Abstract 

We study the field theoretical model of a real scalar field in presence of spacial 
D I inhomogeneity in form of a finite width mirror (material layer). The interaction of 

7^ , ■ the scalar field with the defect is described with position-dependent mass term. We 

calculate the propagator of the theory, the Casimir energy and the pressure on the 
boundaries of the layer. We discuss the renormalization procedure for the model in 
■ dimensional regularization. 

o ■ 

^ ■ 1 Introduction 

00 I Quantum Field Theory (QFT) was developed in the middle of the last century as a theory 

O ' of interaction of elementary particles in otherwise empty, homogenous infinite space-time 

[1]. On the other hand, from the very beginning it was clear that presence of boundaries, 
^ . non-zero curvature or nontrivial topology of the space-time manyfold should influence the 

^ ! spectrum and dynamics of the excited states of the model as well as the properties of the 

ground state (vacuum). 

The first quantitative description of such changes in the vacuum properties was made 
by H. Casimir in 1948. He predicted [2] macroscopical attractive force between two 
uncharged conducting plates placed in vacuum. The force appears due to the influence of 
the boundary conditions on the electromagnetic quantum vacuum fluctuations. Nowadays 
the Casimir effect is verified by experiments with the precision of 0.5% (see [19] for a 
review) . 

The properties of the vacuum fluctuations in curved spaces, investigation of scalar 
field models with various boundary conditions and their application to the description of 
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real electromagnetic effects were actively studied through the last decades, see discussion 
and references in [19], [20] . 

However, it was well understood that boundary conditions must be considered just as 
an approximate description of complex interaction of quantum fields with the matter. A 
generalization of the boundary conditions method has been proposed by Symanzik |3j|. In 
the framework of path integral formalism he showed that presence of material boundaries 
(two dimensional defects) in the system can be modeled with a surface term added to 
the action functional. Such singular potentials with 5-function profile concentrated on 
the defect surface reproduce some simple boundary conditions (namely Dirichlet and 
Neumann ones) in the strong coupling limit. The additional action of the defect should 
not violate basic principles of the bulk model such as gauge invariance (if applicable), 
locality and renormalizability. 

The QFT systems with 5-potentials are mostly investigated for scalar fields. In [1]- 
[6] the Symanzik approach was for the first time used to describe similar problems in 
complete quantum electrodynamics (QED), and all (5-potentials consistent with QED 
basic principles were constructed. 

It seems quite natural to try applying the same method for description of interaction 
of quantum fields with bulk macroscopic inhomogeneities (slabs, finite width mirrors, etc) 
and to study Casimir effects in system of such a kind. There were different attempts to 
quantize electrodynamics in presence of dielectric media (i.e. volume inhomogeneities of 
special kind) see, for instance, [8], [9], none of them was truly successful. The Symanzik's 
method was used to model the interaction of quantum fields with bulk defects in a number 
of papers (e.g. [I0],[T2]-[T5], and others). However most of them were devoted to study 
of a limiting procedure of transition from a bulk potential of the defect to the surface 
^-potential as in [12] . On the other hand, results for the Casimir energy of a single planar 
layer of finite width i are contradictory. Thus, the formulae presented recently in [TB] 
does not coincide with previous calculations made in lilOj. Moreover, the only attempt to 
calculate the propagator in such system was undertaken in [17] where hardly any explicit 
formulae were after all presented. 

Thus, one can see that the specificity of finite volume effects generated by inhomo- 
geneities in QFT has not been yet adequately explored. Our work is dedicated to clarify 
the problem, and to solve existing controversy within an accurate and unambiguous ap- 
proach. We consider a model of massive scalar field interacting with volume defect (finite 
width slab), calculate the modified propagator of the field, the Casimir energy of the slab 
and discuss its physical meaning. 

2 Statement of problem 

Let us consider a model of a real scalar field interacting with a volume defect. In the 
simplest case such defect could be considered as homogenous and isotropic infinite plane 
layer of the thickness £, placed in the X1X2 plane. Generalizing the Simanzik approach, 
we describe the interaction of quantum fields with matter by introducing into the action 
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of the model an additional mass term which is non-zero only inside the defect 

S = ^Jd^x {(j)ix){-dl + m^)(l){x) + Xe{£, xs)(l)\x)) (1) 

where = d'^/dx^ + . . . + 9^/9a;|3)- The distribution function 6'(£, X3) is equal to l/£ 
when Ixal < i/2, and is zero otherwise, in terms of the Heaviside step-function we can 
write it as 6'(£, X3) = [^(^3 + £/2) — 6{x^ — £/2)]/£. Such kind of potential is also called 
patchwise (or piecewise) constant one. In the framework of QFT it was considered for 
the first time in [TD], and later in |13j-|16j. 

To describe all physical properties of the systems it is sufficient to calculate the gen- 
erating functional for Green's functions 



G[J] = N J Dcj) exp{-5[0] + J0}, A^-^ ^ J exp{-So[0] + J(f)} 



(2) 



where J is an external source, S'o(0) = S{(f))\x=o, and normalization for the generating 
functional we have chosen in such a way that G'[0]|a=o = 1- 

Introducing in ([2]) auxiliary fields defined in the volume of the defect only, we can 
present the defect contribution to G[J] as 

exp < — — dx dx3(f)'^{x) r ~ ^ / Dip exp < dx dx^ {~~^ + iy/^ip(f) 

' (3) 

where C is an appropriate normalization constant, and k = X/i. 

With help of projector onto the volume of defect O = 9{xs + i/2) — 6{xz — i/2) acting 

as 

tpOcj) = dx dxsipxj), 
J J-e/2 

we can perform functional integration over 0, and consequently over ip. As a result we 
get 

G[J] = [Detg]-i/2gij5J^ S = D- k{DO)Q-\OD), (4) 
Q = 1 + k{ODO). (5) 

Here the unity operator 1, as well as the whole Q, is defined in the volume of the defect 
only {-i/2, i/2) x and D = (-9^ + m^)'^ is the standard propagator of free scalar 
field. We shall note here that the outlook of (jl]) completely coincides with expression for 
generating functional G[J] in the case of delta-potential term instead of patchwise constant 
one. It is also evident that a straightforward generalization is possible for non-constant 
k(= \/i) with A depending on X3. 

In this paper we calculate explicitly both the modified propagator of the system and its 
Casimir energy, and reveal their dependence on the parameter A describing the material 
properties of the homogeneous defect layer and its thickness i. 



^We operate in Euclidian version of the theory which appears to be more convenient for calculations. 
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3 Calculation of the propagator 

To calculate the propagator 5* defined according to (jlj) let us first derive explicit formula 
for the operator W = Q~^. 

For this purpose we first introduce the Fourier transformation of the coordinates par- 
allel to the defect (i.e. xo, xi, X2). Then for the propagator D of the system without a 
defect one can write 



D(x) = / — 4^e 



(27r)3 J (27r) p| + p2 + ' 
integrating over ps with help of the residue theorem we get 



with E = a/p^ + m^. Then we are able to write the defining (operator) equation for W 
as 

W + kVe2W = 1. (7) 

By construction the mixed p-x^ representation of the free scalar propagator Ve^^x^u) = 
Ve2{x — y) is the Green's function of the following ordinary differential operator 

Kv{x,y)=(^-^ + Vyix-y) (8) 

for V = E"^. Multiplying both sides of ([7]) with K^^ and using obvious relation Kv—Kyi = 
V — V we get 

KpU = -K (9) 

where p = ti + E"^ and U = W — 1. 

The general solution to this (inhomogeneous) operator equation can be written as a 
sum of its partial solution and the general solution of its homogeneous version. Then with 
help of Vp one writes for U 



Uix, y) = -KVpix, y) + «(i/)e^v^ + /?(?/)( 



Here a and f3 — arbitrary functions on y. Imposing the symmetry condition U{x, y) = 
U{y, x) we derive that 

U{x, y) = -nVpix, y) + ae^'+^^v^ + h (e^'-^^v^ + e^y-'^^^P) + ce-("+^)^ 

where a, h and c are some constants now. Introducing W = 1 + U into ([Tj) one gets 

U + kVe2{1 + U) = Q. (10) 
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Requiring that this equation is an identity for all x and y (we remind that U = U{x,y)), 
we find for a, 6 and c 



With help of these expressions we can finally derive the explicit formula for the mod- 
ified propagator of the system. From the definitions of 5* and W, and using ffTOl) we can 
write that 

S = il + U)VE2. (12) 

We divide the general expression of 5 = S{p,X2,y3) into four parts according to the 
position of x^, relative to the defect 





Xs < 


-£/2, ys < -i/2 


S-0 1 


X3 < 


-e/2, y,e{-£/2,£/2) 


S-+, 


X-i < 


-e/2, y, > i/2 




xs e 


{-i/2,i/2), y,ei-i/2,i/2) 



s{p,xs,ys) = { r: fo :rz\fo' ' ' as) 



Other cases could be easily derived using the symmetry properties of the propagator. 
Performing necessary integration according to (fT2|) we get 

S- = V^^ixs - y,) + ^^"""'^^^ ~ ^''^) e^(-3+.3) (14) 
= ^e^(-3+^/2)e^VP ((^ _ i^)ev^fe-^/2) + (^ + i^)eVpW2-?/3)) 

Soo = ^ {2k cosh[(x3 + y3)Vp\ + ev^(l^^-^^l-^)(E - ^p)' + ev^(^-l^3-^^l)(E + 

2\/ P 



with ^ defined in ( ITTl) . 

To the best of our knowledge the only attempt to calculate the propagator for such 
system was presented in [1^ where its final expression was given in terms of "coefficients 
of scattering wave functions" of one-dimensional time-dependent Schrodinger equation. 
However, for the explicit formulae for those coefficients the author refers yet to another 
paper [18| (actually, there is also a misprint in the reference number), where the problem 
of electrons scattering in a powerful laser field is considered and corresponding coefficients 
are presented in the from of infinite series of Bessel functions. The result presented in 
(fT^ is in much simpler closed from, and it raises doubts of correctness of calculations 
presented in [17], [T8] . 
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4 The Casimir Energy 



It is well know that the Casimir energy density per unit area of the defect S can be 
presented with the relation 



S 



1 

TS 



\nG[0] 



2TS 



TT\n[Q{x,y)]. 



(15) 



in the second equation we used (jl]). For explicit calculations we first make the Fourier 
transformation as in (EI). Then 



8 = jJL 



2(27r; 



— Trln[g(p; X3,y3)], 



(16) 



where we also introduced dimensional regularization to handle UV- divergencies and an 
auxiliary normalization mass parameter /i. 

Using the definitions of U and Q we can express ^-derivative of the integrand of flTB]) 
in the following form 

dJnQ = VE2W = -^ 



K 



Then for the energy density we get 



S 



4,~d 



TtU. 



(17) 



We have chosen the lower limit of integration over k, to satisfy the energy normalization 
condition S\k.=o = 0. As we show below the integral is convergent at /t = 0. 
The trace of the integral operator U is straightforward 



TtU 



1 1 2 



dx U{x, x) = 2U + 



-£/2 



4asinh(^iyp) — 



;i8) 



where we already used that a = c. Using a and h given in (fTTj) . one easily notes that 
Tr[/ ~ —£k,/{2E) when k — > 0, thus supporting the above statement. 

Next, putting (|T8l) into (|T71) we can compare our result with previous calculations 
performed in [TO], and also recently rederived in [TT]. Instead of explicit /t-integration in 
f fT7|) . we can differentiate the above mentioned result by Bordag with respect to k to see 
immediately that it coincides explicitly with integrand of (fT7|) . Thus, we write for the 
energy 



£ = fi 



A-d 



\d-l 



2(27r 



In 



-£E 



(19) 



To extract the UV divergencies in (i = 4, let's consider those contributions in £ f|T9|) 
that do not converge while integrated over p. We have 



In 



-IE 



{e'^^iE + ^pf - e-'^P{E - ^pf) 



X 

2E 



E^ 



6 



Hence, within dimensional regularization the energy can be represented as follows 



£ £ fin ~l~ 



div 1 



where 



fin 



27r2 



= In 



-2£E 



AE^p 



(20) 



£. 



Xp 



4-d 



div 



d 



d-l 



p 



AE 



1 - 



A 



2(271)^^-1 

The first item £fin is finite and we removed regularization, while £div is divergent but 
trivially depends on the parameters of the theory and auxiliary parameter p. We add now 
to the action of the model a field-independent counter-term 5S of the form 5S = f + gi~^, 
with bare parameters / and g (of mass dimensions two and one correspondingly). It 
allow us to choose these parameters in such way that the renormalized Casimir energy 
£r defined by the full action S + 6S and considered as the function of renormalized 
parameters appears to be finite both in regularized theory, and also after the removing of 
regularization. 

Thus, for the renormalized Casimir energy we obtain the following result 



fin 



(21) 



where finite parameters fr, Qr must be determined with appropriate experiments. 
The Casimir pressure on the slab is then 

_ d£r _ d£fin Qr 



Taking into account the definition of distribution function ^(£,0:3) one can say that the 
derivative is taken here on condition that the amount of matter (effectively described by 
the defect) in the slab is fixed: j dxs6{i, X3) = 1. Alternatively, one can consider the 
density of the matter to be fixed and calculate the pressure under this condition. Then 
the distribution function has a different normalization condition J dx^9{l, X3) = £, which 
is equivalent to the mere change of variables A — > £A in the formula fl20l) . 



5 Conclusion 

We constructed QFT model of the scalar field interacting with the bulk defect concen- 
trated within a slab of finite width ^. The propagator and the vacuum determinant 
(Casimir energy) were calculated. The later one coincides with results obtained in [TO] , 
[Tl| within a different approach, while the explicit formula for the propagator is given 
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for the first time. The Casimir energy is UV divergent and for its regularization we ap- 
phed dimensional regularization. It allowed us to extract the finite part and to construct 
the counter-terms. The renormalization procedure requires generally two normalization 
conditions to fix the values of the counter-terms with the appropriate experiments. It is 
shown that the Casimir pressure in the system can be calculated in two different ways: 
for fixed density of matter and for fixed amount of matter of the slab. 

Similar problems were considered recently in [16j in the framework of massless scalar 
field interacting with a slab (mirror) of general profile. However, the massless limit of our 
result for the Casimir energy of a single slab ( fT9l) differs from one obtained in [16] , Eq. (68) 
for the case of 'piecewise constant' profile (equivalent to our case). As a validity check we 
appeal to the general perturbation theory. Decomposing the generating functional G{J) 
([2]) in a perturbation series in A, one finds that for the massive theory G{0) is analytical 
at A = with irrelevant (geometry independent) linear term. However, it is evident that 
naive perturbation expansion fails for the limit m — 0, alerting us of non-analyticity of 
the vacuum energy at A = 0. Expanding Eq. (68) of [12] in a power series in A one can 
easily see that it is perfectly analytical with non vanishing linear term, and thus does 
not comply with this general argument. At the same time both the massive and massless 
limits of our result f|T71) . which is equivalent to f|T9|) derived independently by two other 
groups, does posses the required (non-) analyticity properties. 

In our work we considered a model of interaction of quantum scalar field with material 
slab assuming A > 0. One must note that with a simple redefinition of the parameters of 
the system under consideration (i.e. A = —2m?t) one can calculate the Casimir energy 
of two semi-infinite slabs separated by a vacuum gap and interacting through a massless 
scalar field. Similar problem in the framework of quantum statistical physics was first 
solved by Lifshitz, [7J. Comparison with Lifshitz formula, and further generalization of 
the method proposed in this paper to the case of QED is the scope of our future work. 
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